Spin Properties of Very Shallow Nitrogen Vacancy Defects in Diamond 
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We investigate spin and optical properties of individual nitrogen-vacancy centers located within 
1-10 nm from the diamond surface. We observe stable defects with a characteristic optically detected 
magnetic resonance spectrum down to lowest depth. We also find a small, but systematic spectral 
broadening for defects shallower than about 2 nm. This broadening is consistent with the presence 
of a surface paramagnetic impurity layer [Tisler et al., ACS Nano 3, 1959 (2009)] largely decoupled 
by motional averaging. The observation of stable and well-behaved defects very close to the surface 
is critical for single-spin sensors and devices requiring nanometer proximity to the target. 



Isolated defect spins in solids, such as phosphorus 
donors in silicon [1] or the nitrogen-vacancy (NV) cen- 
ter in diamond [2], are considered important building 
blocks for future nanoscale devices, governed by quan- 
tum mechanics. In pure materials, defects can be so 
well-decoupled from their solid host that spin states ap- 
proach a stability normally only found in dilute atomic 
vapors, including coherence times of milliseconds to sec- 
onds [3, 4]. Not surprisingly, atomic defects have over 
the last decade attracted increasing attention motivated 
by their potential for spin qubits in quantum informa- 
tion [5, 6] or for ultrasensitive magnetic detectors with 
nanometer spatial resolution [7-9]. 

The central challenge with many of these endeavors is 
to position the defect of interest in close proximity to 
other circuit elements while retaining their well-defined 
properties known from the bulk. On the one hand, close 
proximity is required for strong enough coupling. For 
example, for the direct coupling to nearby spin magnetic 
dipoles — which scales as r -3 , where r is distance — 
efficient coupling is only achieved at nanometer separa- 
tions. Furthermore, for scanning magnetometry applica- 
tions r directly sets the attainable spatial resolution [7]. 
On the other hand, the coupling will almost always hap- 
pen across a material interface and defects will have to 
be located within nanometers from a surface, potentially 
destabilizing the spin and limiting its usefulness. 

Several mechanisms have been found or proposed to af- 
fect the stability of shallow defects. For single donor spins 
in silicon, for example, the nearby Si/Si02 interface was 
shown to decrease spin coherence times even for donors 
tens of nanometer away due to paramagnetic impurities 
present at the interface [10]. Other possible mechanisms 
include electric surface charge or strain fluctuations that 
may disturb defects through Stark and spin-orbit effects, 
or direct ionization [11, 12]. For nitrogen- vacancy cen- 
ters in diamond, on the other hand, rather little is cur- 
rently known about the spin's performance near the sur- 
face. While functional defects have recently been re- 



ported in < 10 nm diameter nanocrystals [13, 14] and 
within 3 — 4nm from bulk diamond surfaces [15] and co- 
herence times T2 of tens of (J8 have been observed for 
defects at ~ 10 nm proximity [16-18], neither a "shal- 
lowest depth" nor the involved destabilizing mechanisms 
are known. Given the fundamental importance of surface 
proximity for applications, it appears imperative to ex- 
perimentally explore the limits to stability of defects at 
very shallow depths. 

Here we report a systematic study of the spin reso- 
nance properties of single NV defect centers down to 
a proximity of about 1 nm. Defects were produced by 
low energy ion implantation (0.4-5 keV) and investigated 
by optically-detected magnetic resonance (ODMR) spec- 
troscopy. We find well-behaved defects exhibiting a nar- 
row electron spin resonance (ESR) spectrum and coher- 
ence times exceeding 10 microseconds down to the shal- 
lowest investigated depths. We also observe extra line 
broadening for defects shallower than 2 nm. This broad- 
ening is compatible with the presence of surface magnetic 
impurities that are mostly decoupled from the NV spin 
by motional averaging. 

A (lOO)-oriented single crystal of ultrapure diamond 
(< 5 ppb N concentration, Element Six) was used as the 
sample for all experiments. One sample face was im- 
planted with 15 N + or 15 N2~ ions at a series of very low 
energies (0.4 — 5keV, in steps of 0.2 keV) and fluences 
(10 10 - 10 14 N/cm 2 ) [19, 20] (see Figure 1). To form NV 
centers, the sample was annealed for 2h at 800 °C and 
p < 2 • 10 -7 mbar. It is expected that nitrogen atoms do 
not diffuse at these temperatures because the activation 
energy is too high [21-23]. The sample was cleaned by 
boiling it for 24h under reflux in a 1:1:1 mixture of sulfu- 
ric, nitric and perchloric acid and thoroughly rinsed with 
purified water [13]. This procedure is known to remove 
any residues (such as graphite) but to leave the diamond 
sp 3 bonding network intact. Acid treatment also leaves 
a well-defined, oxygen-terminated reference surface [13]. 
Additional details on implantation and sample prepara- 
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FIG. 1: (a) Photoluminescence intensity map of the diamond 
surface. Shallowest 0.4-keV defects are visible at the very top. 
Excitation wavelength was 532 nm and detection bandwidth 
630-800 nm for all experiments except Fig. l(c,d). (b) Close 
up of an implantation spot at 5.0 keV (white square in (a)). 
Some single NV centers are encircled, (c) Luminescence in- 
tensity due to NV - (full circles) and NV° (empty circles) as a 
function of ion energy at a fluence of 8 • 10 11 cm -2 , measured 
using two pairs of filters and corrected for spectral overlap 
[12, 23]. (d) Luminescence intensity as a function of ion flu- 
ence at an energy of 5keV. 



tion are given as Supplementary Material [23]. 

We have performed a detailed inspection of the pre- 
pared diamond surface to validate the sample for later 
spin resonance measurements. Surface roughness was 
determined by atomic force microscopy and was found 
to be very low (i rras = 0.38 nm, over a 300 x 300 nm 2 
window) compared to the defect depth ( > 1 nm) . Angle- 
resolved x-ray photo-electron spectroscopy (ARXPS) was 
used to confirm oxygen termination of the surface and the 
absence of significant graphite residue. The absence of 
sp 2 carbon was further corroborated by confocal Raman 
spectroscopy. No difference was found between implanted 
and non-implanted regions. 

A photoluminescence intensity map of the sample is 
shown in Figure 1(a). Most prominently, we observe that 
optically bright NV centers are visible down to the low- 
est implantation energy (0.4 keV). We have estimated 
the depth of these defects using stopping range of ions in 
matter (SRIM) Monte Carlo simulations [23, 24]; for ex- 
ample, an energy of 0.4 keV corresponds to a peak depth 
of 1.1 ± 0.6nm (see scale in Fig.l(c)). The advantage 
of SRIM calculations, which is the primary method to 
determine defect depths in this study, is that they give 
suitable results of ion implantation depths over a wide 
energy range, including very low implantation energies 
[25] . This is important given the current lack of a precise 
experimental method to directly measure surface proxim- 



ity. The drawbacks of SRIM calculations are that the re- 
sults are inherently statistical (which is addressed below 
by collecting statistics on many defects), and important 
biasing effects like channeling are not considered. We 
have analyzed these effects for our study (see Supplemen- 
tary Material) . In particular we find that ion channeling, 
which could lead to depth underestimation by about 2x 
[22, 26], does not occur for 14 N energies below 0.6 keV 
and is only important towards higher energies [23]. The 
lowest energies, where channeling is absent, are the most 
relevant in this study. 

Figures l(c,d) provide additional photoluminescence 
data that further corroborates this picture. We have 
measured the total photoluminescence intensity both as 
a function of energy and ion fluence and determined the 
relative concentrations of NV - and NV centers. In good 
agreement with earlier studies carried out at higher en- 
ergies [19] we find a monotonic decrease in total photo- 
luminescence with decreasing energy. This decrease has 
been attributed to the vacancy-limited formation of NV 
centers [19]. Since the decrease appears to be mostly due 
to a reduction of NV~ one could conjecture the presence 
of a depth threshold below which the negative charge 
state becomes unstable [12]. We have not, however, ob- 
served any photobleaching and we have only seen a few 
rare cases of fluorescence intermittency among the in- 
vestigated single centers [14] that would support such a 
threshold. The presence of a threshold is also incompati- 
ble with the rapid changes in the spin resonance linewidth 
that we see for the lowest energies (see below). 

We now turn to the core part of this Letter which is a 
study and analysis of electron spin resonance (ESR) spec- 
tra as a function of defect depth. ESR measurements are 
carried out using optically detected magnetic resonance 
spectroscopy [27]. For these measurements, the fluores- 
cence intensity from single, isolated NV center is collected 
while slowly sweeping an auxiliary cw microwave field 
across the spin resonance (~ 2.8 GHz) of the electronic 
ground state. Resonant microwaves induce transitions 
between the m s = and m s = +1 (or m s = —1) spin 
sublevels and lead to an up to 30% reduction in fluores- 
cence. We use this feature to map out ESR spectra of 
single defects and to measure their linewidth and coher- 
ence properties. 

Figure 2(a) collects a series of representative ESR spec- 
tra taken on NV centers at different depths. As a key 
feature we observe increased line broadening as NV spins 
are located closer to the surface: For defects deeper than 
1.8 nm (0.8 eV), there is a clear hyperfine splitting due 
to the 15 N nuclear spin, but for the shallower defects at 
1.5 nm (0.6 keV) and 1.1 nm (0.4 keV) the resonances 
become broad and the hyperfine doublet is barely visi- 
ble or entirely unresolvable. This picture of broadened 
lines was consistent among recorded spectra (> 30 in to- 
tal); for example, we did not find any 0.4 keV defects 
with a resolved hyperfine splitting, while most 0.8 keV 
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FIG. 2: Electron spin resonance measurements on shallow 
NV defect spins, (a) Representative spectra (out of > 30 
total) showing increased line broadening at shallow depths. 
Black dots are experimental points and red solid lines are 
Lorentzian fits. Implanted 15 NV (nuclear spin I = 1/2) are 
distinguished from native 14 NV (/ = 1, 99.6% natural abun- 
dance) by the different hyperfine manifold. The magnetic bias 
field is about 20 gauss, (b) ESR linewidth Auj/2tt (half width 
at half height) plotted against surface proximity d. Black dots 
are experimental values obtained from many separately fitted 
curves such as the ones shown in a). Error bars denote stan- 
dard error. The solid red line is a fit to the black dots based 
on Eq. (2) including motional averaging. The dashed red line 
is the static second moment [Eq. (1)] shown for comparison. 
The linewidth of reference 14 NV is also shown, (c) Spectra of 
a 1.1 -nm defect at zero field (left) and at 3.5 gauss bias field 
(right) rule out the presence of significant surface strain or 
charge (see text). 



not, e.g., the sample or experimental parameters. We 
have also measured a few spectra at higher fluence (up 
to 1 • 10 13 cm -1 ) and found that the line broadening did 
not change between implantation spots of the same en- 
ergy but different fluence. Moreover, no increase in 14 NV 
density is seen on or near implanted areas which elimi- 
nates the possibility that 14 NV centers were formed from 
vacancies created during 14 N implantation, in agreement 
with previous reports [22]. Finally, we did not observe 
a line splitting at zero magnetic bias field [Fig. 2(c)], a 
signature for the presence of electric fields [28], thereby 
excluding the presence of significant surface strain or 
charge. 

In the following we attempt to explain the surface- 
induced line broadening by the presence of paramagnetic 
impurities. Surface impurities have been found at sub- 
stantial density for clean, oxygen-terminated nanodia- 
monds [13]. These nanodiamonds were milled from larger 
crystals and underwent the same surface cleaning proce- 
dure, and are thus expected to exhibit a surface chem- 
istry very similar to our diamond substrate. We can 
model surface impurities by assuming a homogeneous, 
two-dimensional dipolar bath of electron spins (S = 1/2) 
with an areal density of pa ~ 10spins/nm 2 [13], similar 
to the sketch in Figure 3. In the following, we use the 
theory of moments [29], but note that a parallel frame- 
work has been developed for T\ and values in the 
context of paramagnetic impurities in the Si/SiC>2 inter- 
face [10]. The second moment (Aaj 2 ) d ({Auj 2 ) d « Aarij, 
where Alj^ is linewidth) of the NV spin resonance at a 
distance d from the spin bath is found to be 
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and virtually all deeper defects showed a clear hyperfine 
doublet. 

We have quantitatively analyzed the linewidth for the 
shallowest defects (where the pronounced changes are 
seen) by collecting and fitting a number of individual 
spectra and averaging the resultant linewidth parameter 
[Fig. 2(b)]. While we have taken additional spectra at 
other (higher) energies that all show clear hyperfine dou- 
blets, these spectra do not have statistical significance 
and are not included in the figure. Spectra are recorded 
at a single fluence (8- 10 11 cm -2 ) and on defects that lie at 
the perimeter of an implantation dot, where the density 
is low enough to optically isolate individual NV centers 
and residual dipolar broadening by N donor electronic 
spins can be excluded [23] . 

Several control measurements were carried out to en- 
sure that the observed line broadening is indeed a re- 
sult of surface proximity. A number of native 14 NV 
spectra was recorded at each investigated implantation 
spot by focusing slightly into the bulk in order to ver- 
ify that broad lines were a property of the defect, and 



where p = An ■ 10" 7 Vs/(Am), y/(2n) = 2.8 • 10 10 Hz/T 
is the surface electron gyromagnetic ratio, and 9 is the 
angle between the NV axis and surface normal [23]. For 
an NV center near a (100) surface (0 = 54.7°) the angular 
factor in (1) is 3 + 2 cos 2 6 + 3 cos 4 6 = 4. The resonance 
linewidth that can be directly compared to the experi- 
ment is then Aw = \J Awq + (Aw 2 )j, where Acjo is the 
intrinsic linewidth (here Aw /(2ty) = 1.1 ± 0.1 MHz). 

Equation (1) describes the second moment of a quasi- 
static spin bath that does not fluctuate on the time scale 
of the NV dephasing time T 2 * « Aw -1 / 2 . This assump- 
tion is likely invalid for surface spins that are not pro- 
tected by the diamond matrix. For example, spin-lattice 
relaxation times T\ observed for paramagnetic impurities 
in amorphous carbon [30] and sintered detonation nan- 
odiamonds [31] are on the order of nanoseconds. Follow- 
ing the work by Kubo and Tomita [32] we can calculate 
a modified second moment that takes into account fast 
fluctuations ("motional averaging"), 



(Aw 2 )' A = ((Ac 2 ) d r c ) 2 , 



(r e « (A^ 2 ) d /2 ) (2) 



4 



outside spin 



Air 



Interface « 



Diamond 



NVspin 

FIG. 3: Sketch of a nitrogen- vacancy defect spin near a para- 
magnetic surface impurity layer, as described in the text, d 
and denote distance and orientation, respectively, of the de- 
fect spin to the interface plane, and quantization is along the 
NV axis. A potential target spin, representative for sensing 
applications, is also shown. 
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is the correlation time of fluc- 



where t c 
tuations. 

We can compare the model represented by Eqs. (1) 
and (2) through a fit to the experimental data in Figure 
2(b). The three experimental parameters are the sur- 
face spin density pa, the distance d, and the correlation 
time t c . Assuming pA = 10spins/nm 2 [13] and d from 
the SRIM calculation [23] , one obtains a correlation time 
t c < 1 ns. This t c would be rather fast compared to 
the above literature values [30, 31] and typical organic 
radicals. As neither pa nor d are accurately known, our 
value for t c is at best approximate. For example, if we 
allow pa to vary between 0.1 — 10spins/nm 2 and assume 
an underestimation of depth d by up to 2 x , then the 
range of compatible r c varies between 0.1 — 100 ns. Fur- 
ther evidence for a correlation time t c in the nanosecond 
range comes from relaxation time measurements on se- 
lected shallow defects where we observe T\ 3> T\ p « T3 
[33]. 

Finally, we have also measured echo decay times of sev- 
eral 0.4 keV (1.1 nm) defects to establish a lower bound 
for the coherence times T2 of very shallow spins. A rep- 
resentative Hahn echo-decay curve with an echo decay 
time of t = 12 ps is shown in Figure 4; other defects at 
the same depth showed r values between 7 and 12 ps. 
We note that the echo decay profile is Gaussian, which 
is expected if decay is caused by the slowly fluctuating 
13 C nuclear spin environment. If echo decay were set 
by rapidly fluctuating surface impurities, then one would 
rather expect an exponential decay. Consequently, one 
can conclude that surface spins are not the dominant 
source of decoherence and that the limit on T2 implied 
by surface impurities is ~ 10 ps. 

Given the shallow depths of the investigated NV cen- 
ters it is instructive to extrapolate their magnetic mo- 
ment sensitivity, which is the key figure of merit for 
the sensing of external spins and future applications to 
nanoscale magnetic resonance imaging and spectroscopy 
[7]. For example, taking an echo decay time r ~ 12 ps, 
a photon count rate of C = 0.0018 photons/shot and an 
optical contrast between m s = and m s = ±1 states of 
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FIG. 4: Hahn-echo decay of a 1.1-nm 15 NV center and, for 
comparison, for a 7.7-nm 15 NV and a native 14 NV defect. 
Black curves are experimental data and red lines are fits. The 
fit function is given in [23]. 



e = 7% (1.1-nm defect in Fig. 4), we find an optimal ac 
magnetic field sensitivity of -B m ; n ~ (0.57r7e'v / TC) _1 ~ 
2.2//T/\/Hz [34]. Here, the relevant (most suscepti- 
ble) ac frequency is set by the inverse of the echo du- 
ration, i.e., tens of kHz. For a magnetic moment lo- 
cated directly on the surface and taking into account 
the angle of the NV spin, this sensitivity equates to 
a minimum detectable magnetic moment of p m - ln ~ 
Bmin/[(0.96 x p p B ) I '(4vrd 3 )] ~ 0.003 p B /VBz, where p B 
is the Bohr magneton. For dc signal detection, the cor- 
responding magnetic field and moment sensitivities are 
B min « 8Awo/(3V37e v / ^o) ~ 19/xT/vHz and /x min ~ 
0.03 ^b/'v / Hz, respectively, where Iq = 2 • 10 3 photons/s 
is the cw photon count rate and e = 11% [Fig. 2(a)]. 
Even if our depth was underestimated by a factor of 2x 
[23], Pmin would still equate to 0.03/ib/v^Hz (ac) and 
0.2 ^b/V^Hz (dc), respectively. 

In conclusion, we have investigated spin and opti- 
cal properties of single nitrogen vacancy defects in di- 
amond at very shallow depths. Functional defects are 
found down to about 1 nm, and significant broadening 
of the electron spin resonance is only observed for de- 
fects < 2 nm. This surface stability is unmatched by 
other solid-state spin systems, such as phosphorus donors 
in silicon or semiconductor quantum dots, and is a key 
requirement for a number of anticipated quantum and 
sensing applications. In particular, we have inferred a 
sensitivity to outside magnetic moments (such as sur- 
face electron and nuclear spins) that extends down to 
< 0.01 Pb/\/Hz;. To the best of our knowledge this is 
the best magnetic moment sensitivity demonstrated for 
a general-purpose magnetic sensor to date [35]. If com- 
bined with the imaging capabilities of a scanning probe 
apparatus [7, 36], this sensitivity will enable the direct 
mapping of nuclear spins in molecules and thin films with 
chemical specificity and nanometer spatial resolution. 
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^ 1. Ion implantation and annealing 
3 1.1. Ion implantation 
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Ion implantation was carried out with a 5 keV gas-source ion gun (SPECS) combined with a 
Wien mass filter (ExB). Ion beam direction was normal to the diamond surface [a (100) surface] 
with a precision of better than ±1°. Spot size was controlled by a 25 /xm aperture. Fluence was 
measured via the ion current, and adjusted in a way that each spot received the same number of ions. 



< 

Thus, even if spots slightly vary in shape, the number of ions impacting the larger area of the spot 
is the same. Miscalibration of the fluence is the largest error in the photoluminescence intensities 
presented in Fig. 1; in particular, it is most likely responsible for the low photoluminescence of the 
1-keV row. Figure A shows a map of the implantation pattern, energies, and fluences that were 
used on the investigated sample. A detailed description of the ion implantation apparatus is given 
in Ref. [SI]. 

c 

O 1.2. Annealing 

Annealing was carried out in a AJA sputtering machine that allowed heating the diamond sample 
to 800°C while maintaining a pressure < 2 • 10 _7 mbar over the entire annealing duration of 2h, 
followed by a slow overnight cool-down. It is assumed that this pressure is low enough to prevent 
oxidative etching of the top diamond layers. It is also assumed that nitrogen atoms do not diffuse 
under these conditions. Diffusion coefficients for N in diamond have been measured at higher 
temperatures in the context of geology [S2] and were found to follow an Arrhenius relationship, 

O 

^ o o -6.0 eV 

D = 9.7- 10" 8 m 2 /s e k B T . (1) 

> 

At T - 800 °C, the diffusion coefficient is D = 3 ■ 10" 35 m 2 /s. The diffusion length for this D and 
t — 10 4 s is L = \[Dt = 6 • 10 -16 m. This is much smaller than interatomic distances (~ 10 -10 m). It 
is possible that the vacancies created around an implanted N atom will lower the activation energy 
for diffusion, however, such diffusion would be confined to the local area of the N defect and not 
significantly alter the N atoms's position. 

1.3. Exclusion of near-surface 14 NV formation 

There is a possibility that near-surface NV centers may be formed from native 14 N defects due 
to the vacancies created during ion implantation. A signature of such (unwanted) 14 NV formation 
would be an enhanced density of 14 NV centers on or near implantation spots. We have not observed 
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Figure A: (left) Spatial map of ion implantation, (right) Energies and fluences of the various spots. 



such an increase with any of the implantation spots investigated by ODMR. From this we conclude 
that no significant number of near-surface NV centers is formed from native 14 N defects. This is 
plausible given the low intrinsic density of native 14 N of less than < 5ppb [S3]. It is also in agreement 
with earlier studies on N implantation into diamond [S4]. 

2. Estimation of implantation depth for low energy ion implanta- 
tion into diamond 

The depth of defects created by ion implantation depends on the ion species, kinetic energy, 
and the target material. A widely used approach to calculate defect depth is the stopping range of 
ions in matter (SRIM) Monte Carlo simulation [S5]. It has been shown that SRIM simulations give 
suitable results of ion implantation depths over a wide energy range, including very low implantation 
energies [S6]. SRIM is the main approach used to estimate defect depths for our study. Figure B 
shows the calculated defect depth (peak depth) plotted against the ion energy for N ions in the 
range — 5 keV. 
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Figure B: SRIM Monte Carlo simulation for implantation of N ions into diamond. Solid line is peak 
depth (median) and dashed lines are ±1 standard deviation. 



The SRIM simulation does however not consider two potentially important effects that influence 
implantation depth, namely, aperture scattering and ion channeling. These effects are discussed in 
the following. 

2.1. Aperture scattering 

Proximity effects due to scattered ions on the mask are a well known problem in semiconductor 
industry. The scattered ions both reduce the lateral resolution and lead to a depth underestimation 
due to loss of kinetic energy by scattering. The fraction of scattered ions depends on the aperture 
size, material, shape and the ion species and energy [SI]. An upper limit to the fraction of scattered 
ions is given by 
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where R ap is the radius of the aperture and AR the enlargement due to ions that can penetrate the 
aperture and are scattered. For our system R ap = 25 fim and AR f» 5.8 nm (implantation depth of 
5 keV N ions into a Pt aperture curved at 30°) the fraction is / sca tt ~ 10 -3 . For the lowest energies 
used in our experiments the scattering fraction is reduced by an order of magnitude. In addition, 
the aperture is further demagnified by an electrostatic lens to 20 /im and filtered with a second 
200 ixm aperture right in front of the lens. This secondary aperature absorbs ions strongly scattered 
at the main 25 ^m aperature. Thus, we expect no significant impact of aperature scattering on ion 
implantation depth. 

2.2. Channeling 

Ion channeling describes a kinetic ion "traveling" along crystal planes, resulting in a reduced 
number of collisions. Channeling leads to an underestimation of defect depth [S4, S7, S8]. The 
relevant parameters concerning ion channeling are the angle between the ion beam and the crystal 
planes, and ion energy. For a particular crystal orientation, only ions impinging within a certain ac- 
ceptance angle can channel along crystal planes. This angle is called the critical angle. Furthermore, 
channeling does not occur below a certain critical kinetic energy threshold. 
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Figure C: (a) Critical angle and (b) Channeling fraction as a function of ion energy for N implan- 
tation into diamond. The critical energy is 0.6 keV for implantation into the (100) surface. 



We have calculated the critical angle and critical energy for N implantation into (lOO)-oriented 
diamond based on the model by Lindhart [S9, S10]. The results of these calculations are plotted in 
Figure C. Fig. C(a) shows that the critical angle becomes large at low energies and channeling is 
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generally expected under these conditions (the incident angle in our experiment is ±1° off (100)). 
Fig. C(b), on the other hand, shows that the fraction of N ions undergoing channeling is reduced 
at low energies and is zero for ions below a critical energy. This critical energy is 0.6 keV for our 
case. Thus, channeling will occur for a fraction of the ions implanted at higher energy but not for 
the lowest energies (namely 0.4 and 0.6 keV) used in our study. 

The enhancement of penetration depth for the fraction of ions undergoing channeling can only be 
estimated for the low kinetic energies used in our study. Secondary ion mass spectrometry (SIMS) 
measurements on 14 N ions implanted into (100) diamond at higher energies (10-30 keV) showed 
that the deep tail is roughly twice as deep as estimated by SRIM calculations [S4] . Complementary 
molecular dynamics simulations on 4 keV 14 N ions predicted a depth of 14 nm (vs. 6.5 nm by SRIM) 
for channeling ions, underlining the roughly two- fold depth enhancement [Sll]. 

2.3. Cluster implantation 

Cluster implantation is a widely used method in semiconductor industry to increase the through- 
put for shallow implant dopants. The molecule immediately dissociates when hitting the target 
surface and the energy of the atoms is given by the mass ratio of the accelerated molecule [S12]. 
For N2 molecules, it is simply one half the implantation energy. The implantation depth is compa- 
rable to monomer ion implantation [S12]. However, it is known that cluster implantation reduces 
the channeling tail but increases the damage [S12]. Both effects can affect the yield and thus the 
photoluminescence signal; this might be related to the slight step in fluorescence we observe between 
2.6 keV and 2.2 keV implantation spots [Fig. 1(a) in manuscript]. 

3. Photoluminescence measurements 

Photoluminescence and ODMR measurements were carried out on the same home-built inverted 
confocal microscope, except for curves presented in Fig. l(c,d), which used a standard fluorescence 
microscope. NV defects were excited at 532 nm and emitted photons were filtered at an effective 
bandwidth of 630-800 nm and collected by an avalanche photodiode. Single center emission was 
confirmed by photon autocorrelation measurements, ODMR measurements, and the fluorescence 
intensity. A 40x, NA=0.95 air objective (Olympus) was used to focus the laser on to diamond 
sample and collect photons. The sample was mounted to a motorized three-axis stage (Newport 
M-462-XYZ-SD) to navigate over the entire 2x2 surface of crystal. A stationary glass coverslip 
carrying a thin wire was inserted between objective and diamond sample for microwave excitation. 

For selective measurements of NV° and NV - emission two different sets of filters with bandwidths 
of 582-636 nm and 660-735 nm, respectively, were used (Fig. D). Photoluminescence spectra were 
also recorded for selected implantation spots to corroborate the findings from the filter measurements 
(Fig. E). The particular choice of filters is not entirely selective to the two charge states, and the 
curves shown in Fig. l(c,d) have been corrected for the overlap of these filters with the excitation 
spectrum. The raw curves are shown in Fig. D. Corrected intensities were calculated using the 
following equation, 



where the fraction of photons passing the 582-636 nm filter is 0.44 and 0.11 for NV° and NV", 
respectively, and the fraction of photons passing the 660-735 nm filter is 0.26 and 0.66 for NV° and 
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NV , respectively. Numbers are calculated from the spectra given in Ref. [S13]. 
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Figure D: Photoluminescence intensity versus energy and fluence, respectively, measured using two 
sets of filters. 



Photoluminescence versus energy curves that do not discriminate between NV and NV - were 
also measured using the confocal setup for the fluences numbers of 8 • f 11 cm -2 and 1 • 10 13 cm -2 
and gave similar results. 

For ODMR linewidth measurements, laser intensity and microwave power was reduced to ~ 1 /iW 
and a few f 00 kHz Rabi frequency, respectively, to achieve linewidths of < 1 MHz. This linewidth is 
still slightly larger than those reported for intrinsically (nuclear 13 C spin bath) limited high-purity 
samples and it is likely that some residual power broadening is still present. The 1 MHz linewidth, 
however, is sufficient to clearly resolve the hyperfme manifold due to the nitrogen nuclear spin. 



4. Line widths of ODMR spectra 

4.1. Fitting of spectra 

In order to determine an accurate mean number for the linewidth at various implantation energies, 
over 30 spectra were collected and individually fitted. The few spectra that showed a clear, extra 
hyperfme splitting due to proximal 13 C nuclei were excluded from statistics. In the vicinity of each 
implantation spot, 15 NV and 14 NV spectra were deliberatly taken to ensure that line broadening 
effects were indeed characteristic to the spot, and not to local variations in the diamond subtrate. 
We also measured 15 NV spectra at 0.4 keV for different fluences, and found good general agreement 
(unresolvable lines) in all of them. This indicates that there is no qualitative difference between 
ions implanted at different fluence and that there is no dipolar contribution to the linewidth from 
N donor electron spins. 

Line width parameters Aw (0.5 x the full width at half height) for the ODMR spectra were 
estimated by fitting a Lorentzian to each of the hyperfme split lines. For 15 NV , the fit function 
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Figure E: Photoluminescence spectra of selected defects at fluence 1 • 10 13 cm 2 . Ion energy is given 
with each curve. Spectra are vertically offset for clarity. 
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Here, u> is the microwave frequency, and Iq, c, xq, and Aw are free fit parameters that correspond to 
the intensity, the contrast, the center position of the ODMR line, and the linewidth (all in units of 
angular frequency). a/(2ir) = 3.0 MHz is the hyperhne splitting and was assumed a fixed parameter. 

4.2. Calculation of the second moment close to a two-dimensional electronic 
spin bath 

The second moment of the resonance of an electron spin S due to a quasi-continuous, two- 
dimensional layer of electron spins 1—1/2 with density pa (units of spins/m 2 ) is given by 
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(6) 



where 7/ = 7s = 7 = 2tt- 2.8- 10 10 Hz/T is the electron gyromagnetic ratio, fio — in-10~ 7 (Vs)/(Am) 
and r and 9 the distance and angle, respectively, between the NV axis and a surface spin. This 
second moment has units of angular frequency squared. In the following we assume that a small 
bias field Bq is applied along the NV axis and that all surface spins are aligned with this field. 
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For integration we note that 



r(x,y,z) = {x 2 + y 2 + z 2 )i, (7) 

r x z i 

cosfrc, y, z) — — sin(#o) H — cos(#q) (8) 
Lr r J 

where z — d is the distance between spin S and the surface layer, and 9q is the angle between NV 
axis and surface normal. Integration of (6) gives 

, A 2 X l/ig7 4 ^ 2 37rp A [3 + 2cos 2 (e ) + 3cos 4 (^o)] , Q . 

3y 2 p Al 4 h 2 [3 + 2cos 2 (fl ) + 3cos 4 (fl )] 
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For our sample, which has a (100) oriented surface, all NV spin will have the same angle set by 
cos#o = l/v3, or 6*o = 54.7°. For this angle, the second moment becomes 

4.3. Exclusion of line broadening by 15 N spin bath 

We have also calculated the average distance between 15 N defects in the implanted diamond in 
order to exclude extra line broadening due to the N donor electron spin bath. ODMR spectra were 
recorded at a single fluence of 8 • 10 11 cm 2 and on defects at the perimeter of the implantation spots, 
where the density is reduced by about < 1/100 as estimated by fluorescence intensity. Peripheral 
defects were chosen for all measurements in order to resolve individual luminescent centers and to 
ensure N ions were far enough apart. From the ion density we calculate the average (median) N-N 
distance to d ~ 0.5\/0.018 • 10 11 cm 2 « 50 nm. The dipolar interaction at this distance is about 

"D«^i^r«27r.l00Hz, (12) 
4-7T d A 4 

where 7 is the electron gyromagnetic ratio. This value is much less than the base linewidth of 
~ 1 MHz relevant to our study. 

5. Spin echo measurements and fitting of decay times 

Coherence times T2 were measured using a Hahn echo sequence with equal free evolution times 
t' — t before and after the central tt pulse. This pulse sequence produces the decay of the echo 
maximum. The echo decay curves were then fitted by the following equation (adapted from Ref. 
[S14]), 



J(t) = A exp 
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where r c describes the fast initial decay caused by the fluctuating 13 C nuclear spin bath, T2 the 
slower decay of the "echo revivals" appearing at the periodicity r r of the 13 C Larmor precession, 
and n\ = 3, — 2. Here, A, C, r c , T2 and r r are free fit parameters, and N was adjusted to match 
the number of revivals seen. We also tried fitting with arbitrary exponentials m, ri2, but fits would 
either not converge or yielded the same values for r c and T2 within experimental error. Thus, no 
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conclusion can be drawn on the exponent. In principle, the T2 decay exponential will change from 
ri\ — 3 to Tlx — 1 for a rapidly fluctuating environment (such as caused by fast reorientation of 
surface spins), but from the present data we cannot favor one over the other. 
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